On the Use of Transverse Momentum Spectra to probe Thermalisation 
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Abstract 

Transverse momentum spectra of the reaction products from proton-proton and nucleus-nucleus collisions are discussed within 
the phenomenological frameworks of statistical phase-space and string fragmentation models. It will be shown that a string 
fragmentation model can accommodate all of the observations discussed here and consequently it is claimed that no conclusion 
about the thermal nature or collective expansion of the created system may be derived from observed transverse momentum 
spectra. 
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1 Introduction 

With the highly-relativistic heavy-ion physics pro- 
grammes reaching maturity, one can now start to address 
the question of what are the observables to provide un- 
ambiguous evidence for the formation of a Quark-Gluon 
Plasma (QGP) state. 

Before answering this question, physicists have first tried 
to see whether the current data indicate the formation 
of a thermalised system, the logical implication of QGP 
formation. In doing so, transverse momenta of the par- 
ticles produced in the collisions have been investigated. 

In very energetic heavy-ion collisions a large number 
of particles is produced in a relatively small volume and 
consequently many secondary collisions between the re- 
action products take place (rescattering). This makes ex- 
periments with highly relativistic heavy-ion beams well 
suited to search for the formation of a thermally equili- 
brated system. 

Several attemps [1,2] have been made to investigate 
whether the spectra and (relative) yields of the produced 
particles indicate the formation of a thermally and/or 
chemically equilibrated system. 

In addition, interpretation of the observed spectra within 
the framework of hydrodynamical evolution of the cre- 
ated system followed by chemical and thermal freeze-out 
of the hadron source has led to the idea of collective ex- 
pansion effects reflected by the produced particles [3,4]. 



In this paper it will be shown that particle produc- 
tion through the mechanism of string fragmentation also 
leads to particle spectra compatible with the observa- 
tions discussed here, without the necessity of collective 
effects. 



2 The thermal model 

A theoretical description of the momentum spectra of 
particles emerging from a thermalised system is obtained 
by assuming that all the constituents obey a Boltzmann 
distribution. 

In the case of formation of such a thermalised system by 
energetic collisions of heavy ions, it is convenient to in- 
troduce the perpendicular and parallel momentum com- 
ponents with respect to the direction of the incoming 
projectile. These observables are indicated by p± and pu , 
respectively, whereas we also introduce the transverse 
mass m± = \/p\ + m 2 . The distribution of the emerg- 
ing particles as a function of these variables is given by 



dN(p, m,T) = A e 



mj_ dmxdpiidy , (1) 



where T represents the temperature of the system and 
A contains the normalisation to the total number N of 
particles as well as the (mass) terms which account for 
the proper dimensions. 
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Using the azimuthal symmetry of the situation de- 
scribed above, integration over p\\ and tp results in the 
following expression for the differential particle yields 



dN( mjL ,T) 
m^dm± 



where K\ represents the modified Besscl function. 



(2) 



Making use of the relations KAx 3> 1) ^ W — e x and 

V 2x 

Ki(x <1)~ x^ 1 we see from Eq. (2) that at high mo- 
menta an essentially exponential behaviour of the parti- 
cle distribution is obtained, whereas at low momenta the 
differential yield saturates at a level proportional to the 
temperature of the system. 

The theoretical particle spectrum for a thermalised sys- 
tem at a temperature of 200 MeV is shown in Fig. 1. 



Thermal distribution 




Fig. 1. Calculated particle yield for a thermal source of 
200 MeV temperature according to Eq. (2). 

According to the theoretical description of the particle 
spectra produced by a thermalised system, the following 
observations have been made : 

• Various particle species will show the same shape in 
the transverse mass spectra, the so-called m± scaling. 

• At high p± initial hard scattering of the constituents 
starts to play a role [5] and as such an enhancement at 
high p± is expected, the so-called Cronin effect. 

• At low p± one starts to become sensitive to the prod- 
ucts of resonance decays, which induces a rapidity- 
dependent enhancement [6]. 

From this it is seen that, in comparison with experimental 
results, the description of the particle spectra as outlined 
above is only valid in a limited m± interval. 

Furthermore, in general Eq. (2) can not be used in 
comparison with experimental data because of the fact 
that, due to the limited phase-space coverage of the ex- 
perimental apparatus, only a limited range in p\\ can be 
addressed. 

To account for this effect, it is more convenient to express 



the particle yields in terms of the particle rapidity 

'E + i 



y = 0.5 In 



E- 



P\\ 



(3) 



where E represents the particle energy. 
Using the identity pn = m± sinh(y) the Boltzmann dis- 
tribution of Eq. (1) can be written as : 

dN(p, m, T) ~ e~ m ± <*» h <.v)/ T m 2 ± cosh ( y ) dm ± dyd(p. (4) 

Here it should be noted that y represents the rapidity 
in the restframe of the thermalised system, since all the 
descriptions so far have been performed in this frame. 

Unfortunately, the particle distribution of Eq. (4) can 
not be evaluated analytically to provide differential par- 
ticle yields as a function of m± and T within a certain 
rapidity range. For each experimental situation numeri- 
cal integration over the corresponding rapidity coverage 
has to be performed in order to compare the model with 
the experimental data. 



3 What do the data tell us ? 

In this section the theoretical model described be- 
fore will be tested against experimental data. For this 
purpose a selection of data obtained from 208 Pb+ 208 Pb 
collisions at the CERN-SPS with a beam energy of 
158A GeV and p+p collisions at the CERN-ISR at </s 
of 23 GeV has been used. 

Here it should be realised that the ^fs per nucleon- 
nucleon collision in case of the CERN-SPS heavy-ion 
collisions is comparable to that of the CERN-ISR 
proton-proton collisions. 

Concerning the CERN-SPS heavy-ion data we will in- 
vestigate the charged hadron yields as observed in the 
NA49 experiment [7] and the n° yields as observed in 
the WA98 experiment [8] . For both experiments the mid- 
rapidity value in the laboratory system was 2.9. 
The various spectra are shown in Figs. 2 and 3, respec- 
tively. The fact that the spectrum of Fig. 3 can also be 
compared directly with Eq. (4) is easily demonstrated by 
using the identities E = m± cosh(y) and dp\\ = Edy. 

Comparison of the observed spectra with our ther- 
mal model calculations reflected by Eq. (4) leads to the 
preliminary conclusion that indeed the spectra exhibit 
a thermal character compatible with a temperature of 
205 MeV. Furthermore, it has been observed [4,7] that 
the spectra of various particle species exhibit larger 
inverse slopes with increasing particle mass. This phe- 
nomenon has been interpreted in terms of a collective 
expansion of the hadronic source [3,4,7]. 
But is it justified, based on these observations, to con- 
clude that the observed spectra are the result of a 
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Fig. 2. Charged hadron yields [7] as observed within the NA49 
experiment. The solid line indicates a fit with the thermal 
model from Eq. (4). 
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Fig. 3. WA98 observed n transverse mass spectrum [8]. The 
solid line indicates a fit with the thermal model from Eq. (4) . 

thermalised system and that the difference in (inverse) 
slopes for the various particle species is a consequence 
of collective behaviour ? 

To obtain an answer to the above questions, the 
charged pion and kaon yields at mid-rapidity as observed 
in proton-proton collisions at the CERN-ISR [9] have 
been investigated as is shown in Fig. 4. 

From the CERN-ISR data it is seen that the pion 
and kaon spectra (similar results were observed for ir~ 
and K~ spectra) produced in proton-proton reactions 
at similar collision energies as the CERN-SPS heavy-ion 
data are also compatible with our calculated spectrum 
reflected by Eq. (4). In addition, also here it is seen that 
particles with larger mass exhibit a larger inverse slope 
in the m± spectrum. The temperatures obtained here are 
190 MeV and 223 MeV for the pions and kaons respec- 
tively, which even quantitatively agree with the observa- 
tions in the case of heavy-ion collisions [4] . It should be 
noted that in [4] a slightly different fit was used, but in 
the m± range considered here this leads to only marginal 
differences in the extracted temperatures. 



jt + and K* production in p-p 
at 23 GeV cms energy 




Fig. 4. Charged pion (•) and kaon (■) yields observed in p-p 
collisions at y/s = 23 GeV at the CERN-ISR [9]. The solid 
lines indicate fits with the thermal model from Eq. (4). 

Furthermore, even for high energy e + e~ annihilation 
reactions at the PETRA and LEP accelerators [10] the 
particle spectra are found to match calculations based on 
a thermal model. Since it is very unlikely that in proton- 
proton and e+e~ collisions a thermalised system will be 
created and collective expansion effects are present, the 
interpretation of the observed m± spectra in terms of 
a thermalised and expanding hadronic source becomes 
doubtful. 

Further investigations to obtain a more detailed under- 
standing of the underlying reaction dynamics are there- 
fore needed. Below, an attempt is made to explain the 
observed distributions based on a common underlying 
production mechanism. 

4 How to understand the observations ? 

Due to the fact that the interactions described here 
are not calculable within QCD, the interpretation of ex- 
perimental results has to rely on phenomenological ap- 
proaches. 

4-1 Statistical approach 

First attempts to interpret the experimental data of 
high energy nuclear collisions were performed by trying 
to describe statistically the stage in which excitation of 
the incoming matter takes place. However, in their orig- 
inal formulation [11] these descriptions of this so-called 
'early stage' of the interaction failed to reproduce the 
experimental results. 

With the advent of more and more high quality data 
it turned out that after some generalisations had been 
applied to the original ideas, a statistical approach of 
the early stage of the interaction [12] gives rather good 
agreement with the experimental results. 
Furthermore, recently it has been demonstrated [1,13,14] 
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that also the hadronisation process, which follows after 
the expansion and cooling of the early stage, can be de- 
scribed successfully by a statistical model. 

Provided an extrapolation to a full phase-space cover- 
age is performed, these statistical models describe suc- 
cessfully the various hadron abundancies (i.e. the chem- 
ical composition) of the interactions as observed in the 
various experiments. However, they provide no descrip- 
tion of the dynamical properties of the created system as 
for instance reflected in the m± spectra of the reaction 
products. 

These models also fail to describe the observed event 
topologies like the leading particle effect and jet struc- 
tures in e + e~ annihilations and in addition need the ex- 
plicit assumption that the temperature of the system re- 
mains constant during the complete hadronisation pro- 
ces. Furthermore, an ad-hoc strangeness suppression fac- 
tor (7 S ) has to be introduced to accommodate the ob- 
served yields of particles containing strangeness. 

All this makes one wonder whether such statistical ap- 
proaches are the right way to go. Perhaps a description of 
the interaction process on a much more elementary level 
could provide a proper description of all the observed 
characteristics. 



4-2 String fragmentation model 

Based on an intuitive interpretation of the Dual Res- 
onance Model a string model of hadrons has been devel- 
oped [15] . This model was seen to provide a simple expla- 
nation for confinement and also leads to the fact that the 
particle spin is proportional to the square of the particle 
mass, the so-called Regge behaviour. 

Based on the success of the string model, such a string 
formalism was also used to provide a description of high 
energy collisions [16,17] and it was seen that such a phe- 
nomenological approach provided a quantitative descrip- 
tion of experimentally observed effects such as : 

• Leading particle effect and jet structures. 

• The existence of a rapidity plateau around mid- 
rapidity. 

• Limiting fragmentation, resulting in a logarithmic 
growth of the particle multiplicity as a function of the 
collision energy. 

• Relatively small values of the transverse momentum 
of the produced particles. 

Independently of the classical string phenomenology, a 
QCD inspired string fragmentation model was developed 
[18] which also provided a quantitative description of 
high energy interactions. The latter description has the 
advantage that it is well suited to be implemented in com- 
puter simulation programs, the so-called Lund Monte 
Carlo program. 



Both scenarios turn out to give similar results and in this 
paper the Lund formalism will be used. 



4-2.1 Transverse momentum spectra 

In section 5 of [18] it is demonstrated that fragmen- 
tation of a stretched-out colour field provides a common 
mechanism to successfully describe particle production 
for all of the following processes : 

• Quark fragmentation in e + e~ annihilations. 

• Meson fragmentation in hadronic collisions. 

• Proton fragmentation in hadronic collisions. 

• Proton fragmentation in deep inelastic scattering 
(diquark fragmentation) . 

In view of the similarities concerning the various m± 
spectra mentioned before, it is therefore tempting to 
try to explain the observed mi distributions within the 
framework of a string fragmentation model. 

As outlined in [18], hadron-hadron interactions can be 
described by representing the initial hadrons as massless 
quarks connected by a string in "yo-yo mode" for a me- 
son and by a Y-shaped string for a baryon. Due to the 
relative motion of the quarks in the initial hadron the 
colour field will not be completely straight, but effectively 
in the overall cms of the collision a linear description is 
obtained by projection on the direction Q + k, where Q 
represents the momentum transfer and k the primordial 
momentum of the struck (anti)quark. 
Comparison of the string fragmentation model with var- 
ious data indicates that the value of k± (with respect to 
the direction Q + k) is of the order of about 400 MeV/c. 
The remaining momentum — k of the system when the 
reaction takes place is distributed among the particles 
in the " target" fragmentation region in which a particle 
with a large energy-momentum fraction z also takes a 
large fraction of this remaining momentum. 

Limiting ourselves to soft processes and light-particle 
production at mid-rapidity (i.e. low z values), it follows 
that the effects of the primordial k± arc marginal in 
our discussion concerning the observed m± spectra men- 
tioned before. In particular for initial baryons one quark 
will always be close to the junction which leads to an es- 
sentially linear picture describing the diquark fragmen- 
tation along the incoming beam direction. 

In the following the string fragmentation model will 
be used to derive the m± spectra for mesons produced 
at mid-rapidity in the case of hadron fragmentation. 
Consider the case where, after string fragmentation, a 
quark of flavour a and an anti-quark of flavour f3 com- 
bine into a meson with transverse mass m± and energy- 
momentum fraction z. According to [18] the probability 
for this process to occur in the case of symmetric string 
fragmentation is given by the so-called scaling function 
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f a p{z, m ±) which is found to be 

f af3 (z,m ± ) = N afj (l-zr-e- bm ^ z , 



(5) 



where a and b are free parameters and N a p represents a 
normalisation constant. 

For low z values the probability reflected in Eq. (5) can 
be written as 



f a p{z,m±) 



-bm\l z 



(6) 



whereas for z — > f the expression of Eq. (5) basically 
becomes (1 — z) a . 

Following the phenomenology outlined in detail in sec- 
tion 2.6 of [18] it is seen that the parameters a and b are 
strongly correlated, since increasing a or decreasing b will 
increase the mean multiplicity and affect the shape of 
the transverse momentum spectra. In order to match the 
experimentally observed spectra and multiplicities only 
values along a curve in the (a, b) plane are allowed. 
Good agreement with various data is obtained by a = 1 
and b — 0.4 GeV~ 2 . Using these values of a and b, it 
is seen that particle production at high z values is sup- 
pressed and that the distribution at low z values peaks at 
z = 0.4m5_ (for low mj_ values) and vanishes for smaller 
z values. 

As discussed before, light-particle production at mid- 
rapidity due to soft processes involves low z values. If 
we thus limit ourselves to this part of phase-space and 
make use of E = mj_cosh(y), p« = m± sinh(y) and 



= (E±P\\ 



i hadron 



/(E ±p\\) total, the latter of which 



are the invariant lightcone variables describing the "for- 
ward" and "backward" fragmentation in case of finite jet 
energies [18], it is readily seen that in our case z ~ mj_. 

This eventually leads to the following distribution 
function concerning light-particle production in the 
central region : 



1 

m± 



(7) 



where A comprises both the constant b and the term 
(E ±p\\)total corresponding to the forward or backward 
fragmentation. 

From this probability the distribution of the particles 
emerging around mid-rapidity is found to be 



dN(p, m) <~ e Ami mi cosh(j/) &m±dydip. 



(8) 



Here it is seen that also the phenomenological string frag- 
mentation model results in a transverse mass distribution 
of the particles produced in high energy collisions which 
exhibits a shape compatible with a thermal source. 



Obviously one could now raise the question of how Eq. (8) 
would fit the experimental spectra. 

Taking again the charged pion and kaon yields from 
the CERN-ISR proton-proton data as shown in Fig. 4 
and performing a fit of the data with Eq. (8) shows that 
the model describes the data rather well as can be seen 
from Fig. 5. 



it + and K* production in p-p 
at 23 GeV cms energy 




Fig. 5. Charged pion (•) and kaon (■) yields observed in p-p 
collisions at y/s = 23 GeV at the CERN-ISR [9] fitted against 
the string fragmentation model from Eq. (8) . 

Similar results are obtained for the tt~ and K~ data. 
The inverse slopes of the fit, which in a thermal model 
would represent the temperature of the system, are found 
to be 218 MeV and 257 MeV for the pions and kaons, 
respectively. As outlined in section 2.6 of [18], also from 
the string fragmentation model it is seen that the inverse 
slope becomes larger for heavier particles since they need 
a slightly higher energy-momentum fraction z at their 
creation. Considering momentum values in the range of 
0.5-1 GcV/c it is seen that kaons exhibit an increase in z 
of about 20% compared to pions, resulting in a compa- 
rable increase in inverse slope of the m± spectra. This is 
perfectly in line with the observations in the case of the 
ISR data discussed above and also in the case of heavy- 
ion collisions at the CERN-SPS [4]. 
Here it should be noted that the production of particles 
containing strangeness needs a separate study concern- 
ing heavy-quark production in the string fragmentation 
model. Since this is outside the scope of this paper we 
will not elaborate on this subject here, but it is seen [18] 
that, apart from a relative suppression factor of about 
■| for strange-quark production, the fragmentation into 
kaons proceeds similarly to that into pions. 
It is thus seen that the string model also automatically 
yields a suppression of heavy-quark production which 
turns out to be in agreement with experimental observa- 
tions and even matches the ad-hoc strangeness suppres- 
sion factor 7„ = 0.3 which had to be introduced in the 
statistical model [1,13,14]. 
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Comparison of the nucleus- nucleus data from Figs. 2 
and 3 with the string fragmentation model of Eq. (8) is 
presented in the Figs. 6 and 7 below. 
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Fig. 6. Charged hadron yields [7] as observed within the 
NA49 experiment. The solid line indicates a fit with the string 
fragmentation model from Eq. (8). 




: , i , i , i , i , i , i ± L 

0.5 1 1.5 2 2.5 3 3.5 

m ± (GeV) 

Fig. 7. WA98 observed 7r° transverse mass spectrum [8]. The 
solid line indicates a fit with the string fragmentation model 
from Eq. (8). 

Also in the case of the m± spectra of the hadrons pro- 
duced in nucleus-nucleus collisions the string fragmenta- 
tion model is seen to describe the data rather well. The 
inverse slopes from the NA49 and WA98 data are found 
to be 237 MeV and 227 MeV respectively. 
In order to compare the various observed (inverse) slopes, 
the values obtained from the NA49 and WA98 data have 
to be corrected for the rapidity coverage. Using the values 
of the specified rapidity intervals expressed in the overall 
cms of the collision and following the arguments given 
above, both inverse slopes amount to about 200 MeV af- 
ter the rapidity correction. 

From this it is seen that indeed the string fragmentation 
model is able to describe all the data discussed in this 
paper in a consistent way. 



5 Conclusions 

Various models describing particle production were 
confronted with data from proton-proton and nucleus- 
nucleus data. 

It was seen that both a thermal model and a phenomeno- 
logical string fragmentation model describe the observed 
transverse momentum spectra, where the string fragmen- 
tation model automatically comprises the increase in in- 
verse slope with particle mass as has been observed in 
both the p-p and heavy-ion data discussed here. To ac- 
commodate this in the thermal model, different freeze- 
out temperatures have to be assumed for the various par- 
ticle species or a collective expansion of the system has to 
be imposed, in which the expansion velocity is the same 
for all particle species. 

Concerning the production of various quark flavours, 
the string fragmentation model was seen to yield auto- 
matically suppression of heavy quarks as compared to 
the light u and d flavours. The strangeness suppression 
factor of about i turned out to match the ad-hoc intro- 
duced 7 S from statistical models. 

Finally it should be noted that the string model is the 
only one which accounts for the observed event topolo- 
gies like for instance the leading particle effect and jet 
structures. 

From this it is concluded that a phenomenological 
string fragmentation model describes all aspects of the 
data discussed here and consequently it is not justified 
to draw any conclusions concerning a thermal nature or 
a collective expansion of the created system from the ob- 
served transverse momentum spectra. 
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